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(S) Defective pixel correction circuit 



(5?) A defective pixel correction circuit corrects a 
defective pixel in a solid imaging device such as 
a CCD exactiy and sufficiently. A boundary 
detection circuit (11) calculates magnitudes of 
boundaries from . signals of eight peripheral 
pixels taken in a pixel taking-in circuit (10) and a 
boundary ordering circuit (12) compares the 
calculated magnitudes of the boundaries with 
one another to order the magnitudes of the 
boundaries. An interpolation circuit (14) pro- 
duces an interpolation signal in accordance 
with an inter-polation method determined by an 
interpolation method determining circuit (13) 
on the basis of the order of the ordered bound- 
aries to correct the defective pixel . 
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BACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION 

5 The present invention relates to a defective pixel correction circuit for correcting a defective pixel in a solid 

imaging device such as a CCD. 

DESCRIPTION OF THE PRIOR ART 

w Correction of a defective pixel in a CCD or the like is heretofore performed by a so-called 0-order hold 

interpolation method in which positional information of a defective pixel is previously stored and the defective 
pixel is simply replaced by information obtained from a pixel adjacent to the defective pixel on the basis of the 
positional information. In this correction method, however, as shown in Fig. 1, when a boundary between light 
and dark portions of an image is present just before a defective pixel, for example, the place X in which the 

is defective pixel is corrected is conspicuous. 

In order to solve such a problem, JP-A-3-296374 proposes a correction method in which levels of image 
signals produced from a plurality of pixels in the vicinity of a defective pixel are compared with one another to 
select image information for the pixel used for interpolation on the basis of the comparison result so that image 
information at the position of the defective pixel is interpolated. The patent publication as mentioned above 

20 discloses the following measures. 

More particularly, as shown in Fig. 2, absolute values of signal level differences between adjacent pixels 
(boundaries I, II, III and IV) for five pixels including pixels D and E disposed before and after and adjacent to 
a defective pixel X and pixels A, B and C disposed before one horizontal line are calculated. 

The absolute values of the output signal level differences between the adjacent pixels are compared be- 

25 tween the boundaries I and ll t I and 111, IV and II, and IV and III to produce four kinds of control signals a to d 
determined in accordance with Table 1 on the basis of the comparison among the four boundaries. 
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As shown in Figs. 3Aand 3B, the control signals a to d correspond to four kinds of stripe patterns including 
a horizontal pattern shown in Fig. 3A, a right tilted pattern shown in Fig. 3B, a left tilted pattern shown in Fig. 
2C and a vertical pattern shown in Fig. 3D. The d fective pixel is replaced by an average of the output data 
55 of th pixels D and E for the horizontal pattern of Fig. 3 A, the output data of the pixel E for th right tilted pattern 
of Fig. 3B, th output data of the pixel D for the left tilted pattern of Fig. 3C and the output data of the pixel B 
for the vertical pattern of Fig. 3D to b corrected. 

In the correction circuit proposed in the above-mentioned JP-A-3-296374, how ver, the absolute values 
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of the level differences between the pixels are calculated from the levels of the image signals produced from 
the f iv pixels including the right and left pixels of the defective pixel and the three pixels disposed b fore one 
line of the defective pixel and are compared with one another to thereby identify a pattern of an object to be 
imaged in the vicinity of the def ctive pixel and accordingly the pattern is sometimes identified in error. 
5 For example, in Fig. 4A, when outputs of peripheral pixels ar A=100, B=120\ C=200, D=110 and E=215, 

absolute values of boundaries are calculated as follows: 

Boundary I = | A- D| = 10 
Boundary II = | B - A | = 20 
Boundary III = | C - B I = 80 
10 Boundary IV = |E-C| =15 

Accordingly, this case corresponds to No. 14 of Table 1 and the control signal d is selected. However, since 
the pattern of Fig. 4A is tilted right, the pattern must be identified as the right tilted pattern of Fig. 3B and the 
control signal b must be selected. Further, in Fig. 4B, when outputs of peripheral pixels are A=100, B=110, 
C=200, D=120 and E=205, absolute values of boundaries are calculated as follows: 
15 Boundary I = I A - D | =20 

Boundary II = I B - A I = 10 
Boundary III = |C- B| =90 
Boundary IV = |E-C| =5 
Accordingly, this case corresponds to No. 7 of Table 1 and the control signal a is selected. However, since the 
20 pattern of Fig. 4B is vertical, the pattern must be identified as the vertical pattern of Fig. 3D and the control 
signal d must be selected. 

As apparent from the above two examples, there is a problem that the correction circuit cannot identify 
the pattern at the boundary in the vicinity of the defective pixel exactly. Further, since the patterns of an object 
to be identified are limited to four kinds of patterns including the horizontal, left tilted, right tilted and vertical 
25 patterns, sufficient and fine correction is impossible. 

SUMMARY OF THE INVENTION 

It is an object of the present invention to solve the above problems in the prior art by providing an excellent 

30 defective pixel correction circuit capable of correcting a defective pixel exactly even when a complicated bound- 
ary is present in an object to be imaged in the-vicinity of the. defective pixel. 

In order to achieve the above object, the defective pixel correction circuit according to the present invention 
comprises a pixel taking-in circuit for taking in signals produced from pixels in the vicinity of a defective pixel 
supplied from a digital image, a boundary detection circuit for calculating magnitudes of boundaries between 

35 peripheral adjacent pixels from the signals produced from the taken-in pixels, a boundary ordering circuit for 
comparing the magnitudes of the detected boundaries to order the boundaries, an interpolation method deter- 
mining circuit for determining an interpolation method on the basis of inputs from the boundary ordering circuit, 
and an interpolation circuit for producing an interpolation output in accordance with the interpolation method 
determined by said interpolation method determining circuit. 

40 With the above configuration, according to the present invention, the magnitudes of the boundaries be- 

tween the peripheral pixels are compared with one another to order the magnitudes of the boundaries so that 
even when a complicated boundary is present in the vicinity of the defective pixel, the pattern of the boundary 
can be identified. Further, the interpolation method determining circuit previously determines the optimum in- 
terpolation method for each identified pattern so that correction of the defective pixel can be made exactly for 

45 each boundary pattern. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a diagram for explaining a conventional correction process using a 0-order hold interpolation; 
50 Fig. 2 shows a pixel arrangement in the vicinity of a defective pixel in a prior art; 

Figs. 3A to 3D show stripe patterns identified in a prior art; 
Figs. 4A and 4B show patterns which cannot be identified in a prior art; 

Fig. 5 is a block diagram schematically illustrating an embodiment of a defective pixel correction circuit of 
the present invention; 
55 Fig. 6 shows a window of pixels taken in by a pixel taking-in circuit 

Fig. 7 shows a window of pixels taken in a pixel taking-in circuit in a frame image in the interlaced scanning 
system; 

Fig. 8 shows positions of boundaries; 
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Fig. 9 shows an actual example of d tected boundari s; . . 

Figs. 10Aand 10B show examples of actual boundaries presumed from detected boundari s; 
Figs. 11Aand 11B show examples of interpolation methods suitable for boundary patterns; and 
Figs. 12A to 12L show 12 kinds of interpolation methods. 
5 - 
DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Embodiments of the present invention are now described with reference to the accompanying drawings. 
Fig. 5 is a block diagram schematically illustrating an embodiment of the present invention. Numeral 10 
10 denotes a pixel taking-in circuit for taking in signals of nine pixels corresponding to a window of 3x3 supplied 
from a digital imaging device such as a CCD, numeral 11 denotes a boundary detection circuit for calculating 
magnitudes of boundaries between adjacent pixels on the basis of eight peripheral or surrounding pixels taken 
in by the pixel taking-in circuit 10, numeral 12 denotes a boundary ordering circuit for ordering the boundaries 
on the basis of the magnitudes of the boundaries calculated by the boundary detection circuit 11, numeral 13 
15 denotes an interpolation method determining circuit for determining an interpolation method on the basis of an 
output of the boundary ordering circuit 12, and numeral 14 denotes an interpolation circuit for performing in- 
terpolation using signals from the eight surrounding pixels taken in by the pixel taking-in circuit 10 in accordance 
with the interpolation method determined by the interpolation method determining circuit 13. 

Numeral 15 denotes a memory for storing a position of a defective pixel, numeral 16 denotes, a system 
20 controller for judging whether the position of the center pixel taken in by the pixel taking-in circuit 1 u' is stored 
in the memory 1 5 or not to control a change-over circuit 1 7, and numeral 17 denotes the changing-over circuit 
for selecting between a signal of the center pixel taken in by the pixel taking- in circuit 10 and an interpolation 
output from the interpolation circuit 14 in accordance with a control signal from the system controller 16. 
Operation of the embodiment is now described. 
25 A method of taking in an image is first described. The pixel taking-in circuit 1 0 takes in the signals of nine 

pixels corresponding to the window' of 3x3 and produces the signals. That is, as shown in Fig. 6, the pixel taking- 
in circuit 10 produces a signal x11 from X11 corresponding to the center pixel and signals xOO to x02, x10, x12 
and x20 to x22 from XOO to X02, X1 0, X1 2 and X20 to X22 corresponding to the surrounding or peripheral pixels, 
respectively. However, in the case of an interlaced signal in the interlaced scanning method, interpolation is 
30 made in a field and accordingly the nine pixels correspond to a window of 3x3 in the field. In other words, in 
the case of the interlaced scanning, the pixel taking-in circuit takes in signals from pixels disposed before two 
lines and after two lines in a frame picture as shown in Fig. 7. 

A method of calculating magnitudes of the boundaries between the surrounding pixels is now described. 
The boundary detection circuit 11 calculates absolute values of differences between the signals from adjacent 
35 pixels of the eight surrounding pixels taken in by the pixel taking-in circuit 10 as shown in Fig. 8. However, the 
magnitudes of the boundaries between the vertically adjacent pixels are multiplied by a coefficient k. That is, 
the magnitudes d1 to d8 of the boundaries are calculated by the following equation 1: 

Equation 1: d1 = k- 1x10 - x00| 
d2 = |x00 - x0l| 

40 d3 = 1x01 - x02| 

d4 = k- |x02 - x12| 
d5 = k - 1x12 - x22| 
d6 = Ix22 - x2l| 

d7 = 1x21 - x20| 

45 d8 = k- |x20 - x10| 

where the coefficient k is a positive real number which is used to correct the strength of the correlation between 
pixels in the horizontal and vertical directions caused by different distances between pixels for calculation of 
a signal difference in the horizontal and vertical direction due to the aspect ratio of the image and the interlaced 
scanning. 

so A method of comparing the boundaries in magnitude is now described. The boundary ordering circuit 12 

compares the magnitudes d1 to d8 of the boundaries with one another with respect to all combinations of 56 
sets including d1 and d2, d1 and d3, .... d7 and d8 and produces the order of the boundary magnitude d1 as 
D1, the order of d2 as D2, and the orders of the boundary magnitude d3 to d8 as D3 to D8. respectively. A 
combination of D1 to D8 is named a boundary order pattern. Further, when the comparison results are equal 
55 to each other, the relation in magnitude is determined in accordance with the following rules: 

Rul 1 : when the comparison results in magnitude of the boundaries in the v rtical and horizontal directions 
are equal, it is assumed that the magnitude of the boundary having smaller one of the correlation between 
the horizontally adjacent pix Is and the correlation between the v rtically adjacent pixels is larger than 
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that of the oth r. Further, in the interlaced scanning system, since the correlation betwe n the horizontaliy 
adjacent pixels is generally higher than th correlation betw en the vertically adjacent pixels, the boundary 
in th vertical direction is assumed to be larger in magnitude when the comparison results in magnitude 
of the boundaries in the vertical and horizontal directions are equal. 
s Rule 2: when the comparison results in magnitude of the boundaries in the same direction are equal, va- 

rious rules can be considered since the positional relation is symmetrical, while th boundary having a 
smaller suffix is assumed to be larger than the boundary having a larger suffix in the embodiment. 
The rules 1 and 2 are combined to form the relation of the following equation 2 and accordingly when the 
comparison results are equal, the boundaries are ordered in accordance with the equation 2: 
10 Equation 2: d1 > d4 > d5 > d8 > d2 > d3 > d6 > d7 

The above operation is now described with reference to a definite example. When the relation of d1 to d8 
is defined by the following equation 3, the boundary order pattern is calculated. 

Equation 3: d3 > d5 > d1 > d6 > d8 > d2 = d7 > d4 
That is, d1 is in a third place and accordingly the order D1 of d1 is three. Further, d2 and d7 are both in 
15 the sixth place, while since it is judged that d2 is larger than d7 (d2 > d7) in accordance with the above rule, 
d2 is in the sixth place and d7 is the seventh place. Accordingly, D1 to D8 are determined to be the following 
equation 4: 

Equation 4: D1 = 3, D2 = 6, D3 = 1, D4 = 8, D5 = 2, D6 = 4, D7 = 7, D8 = 5 

A method of determining an interpolation method from the comparison results is now described. The 
20 boundary order pattern produced by the boundary ordering circuit 12 represents a combined order of the mag- 
nitudes of the eight boundaries and accordingly there are eight factorial, that is, 8!=40320 kinds of patterns. 
The interpolation method determining circuit 13 assigns a peculiar number to each of the boundary order pat- 
terns produced by the boundary ordering circuit 12. This number is named the boundary order pattern number. 
The boundary order pattern number is calculated from the boundary order pattern D1 to D8 produced by the 
25 boundary ordering circuit 12 in accordance with the following equation 5. 

Equation 5: PN = D1 + 8(D2 + 8(D3 + 8(D4 + 8(D5 + 8(D6 + 8(D7 + 8 D8)))))) - 

Since D1 to D8 do not have the same value, the pattern number PN does not have a continuous number. 

In order to determine the interpolation method for the boundary order pattern classified above, the inter- 
polation method determining method 13 uniquely sets the interpolation method to be outputted for each of the 
30 patterns previously. The relation of the order of the boundary magnitude in the comparison result and the in- 
terpolation method is now described definitely. 

When it is assumed that the boundary order pattern D1 to D8 produced by the boundary ordering circuit 
12 is set as indicated by the following equation 6 and the order of the boundary magnitude is set as shown in 
Fig. 9, boundaries of an actual peripheral image as shown in Figs. 1 0A and 10B are considered as boundaries 
35 of the actual peripheral image to be presumed from the boundaries shown in Fig. 9. 

Equation 6: D1 = 7, D2 = 8, D3 = 1, D4 = 5, D5 = 4, D6 = 3, D7 = 6, D8 = 2 

Thus, when the patterns of the peripheral image are different, optimum interpolation methods therefor are also 
different. The interpolation method of Fig. 11A is considered to be optimum for the pattern of the peripheral 
image of Fig. 10Aand the interpolation method of Fig. 11B is considered to be optimum for the pattern of Fig. 
40 10B. 

However, when the boundary order pattern as shown in Fig. 9 is obtained by the boundary ordering circuit 
12, it is difficult to uniquely determine the optimum interpolation method from a large number of candidates 
for the interpolation methods including those as shown in Figs. 11 A and 11 B. Accordingly, a simulation is pre- 
viously made in accordance with the following method to determine an optimum interpolation method for each 
45 of the boundary order patterns. 

The simulation uses a digital image having no defective pixel. 

First of all, nine pixels corresponding to a window of 3x3 having a center of coordinates (x,y) are taken in 
from the image. In order to specify the process for the window of 3x3 having the center of coordinates (x,y), 
a value of the center pixel is attached with a suffix of (x.y) such as x11(x,y). Further, in the case of the interlaced 
50 scanning system, pixels of the window of 3x3 for a field image are taken in the same manner as in the pixel 
taking-in circuit 10. 

The nine pixels taken in by the circuit 1 0 are processed by the boundary detection circuit 11 and the bound- 
ary ordering circuit 12 as described above so that the boundary order patterns D1 (x,y) to D8(x,y) are calculated 
and the boundary order patt rn number PN(x,y) is calculated in accordance with the equation 5. 
55 Interpolation outputs interpolated by 12 kinds of int rpolation methods shown in Figs. 12A to 12L are cal- 

culated from th eight surrounding pixels taken in by the circuit 10 so that an interpolation output calculated 
by an interpolation method 1 is assum d to y1(x,y), an int rpolation output by an interpolation method 2 is as- 
sumed to y2(x,y) and interpolation outputs by interpolation methods 3 to 12 are assumed to y3(x,y) to y12(x,y), 
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respectively. Figs. 12Ato 12L show the 12 kinds of simple interpolation methods, while other interpolation meth- 
ods using an average or a weighted average of a plurality of pixels may be considered. Gen rally, an interpo- 
lation output yn can be expressed by the following equation 7: 
Equation 7: yn = alxOO + a2 x01 + a3*x02 + a4-x12 ♦ a5oc22 + a6x21 + a7 x20 + a8x10 

5 where a1 to a8 represent coefficients of actual numbers. 

On the other hand, in this simulation, pixels are taken in from an image having no defective pixel and ac- 
cordingly a signal x11(x,y) produced by a center pixel X11(x,y) in the window of 3x3 can be considered to be 
a target value for correction of the defective pixel. Thus, interpolation errors of the respective interpolation 
methods can be calculated. An interpolation error ERRi(x.y) of an interpolation method i is given by the following 

10 equation 8: 

Equation 8: ERRi(x.y) = I yi(x,y) - x11(x,y) I 

where i = 1, 2, .... 12 

The interpolation error ERRi(x,y) is calculated for the whole image and a squared average of the error is 
calculated for each boundary pattern number. That is, when an interpolation method i is adopted in the bound- 
15 ary order pattern number j, a squared average TOTALERR(ij) of the interpolation error is given by the following 
equation 9: 

Equation 9: TOTALERR(iJ) = £Hj(x,y)ERRi(x,y)ERRi(x,y) (x.y) e all pixels 

where 

Hj(x T y) = 1:(forj=PN(x,y)) 

20 Hj(x,y) = 0: (for j*PN(x,y)) 

By selecting the interpolation method having the minimum squared average of the squared averages TO- 

TALERR(1 j), TOTALERR(2,j) TOTALERR(1 2 j) of the interpolation errors by the 12 kinds of interpolation 

methods, the optimum interpolation method can be determined uniquely for the boundary order pattern number 
j. The above operation is made for all of the boundary order pattern numbers. 

25 Even when it is difficult to uniquely determine the interpolation method suitable for the boundary pattern 

obtained from only the surrounding pixels by means of intuitive judgement, the simulation can be performed 
for image data having various conditions to thereby determine the interpolation method suitable for each bound- 
ary pattern. 

The interpolation method is now described. The interpolation circuit 14 performs calculation for signals pro- 
30 duced from the eight surrounding pixels in accordance with the 12 kinds of interpolation methods (equation 
10 described later) shown in Figs. 12Ato 12Lto produce interpolation signals y1 to y12. An interpolation signal 
of the interpolation method determined by the interpolation method determining circuit 13 is selected from the 
above 12 kinds of interpolation signals and is produced. The interpolation method expressed by the equation 
7 can be added as the interpolation method to be selected except the 12 kinds of interpolation methods shown 
35 in Figs. 12Ato 12L. 

Equation 10: Interpolation Method 1: y1 = x10 
Interpolation Method 2: y2 = x12 
Interpolation Method 3: y3 = x01 
Interpolation Method 4: y4 = x21 
40 Interpolation Method 5: y5 = (x10 + x12)/2 

Interpolation Method 6: y6 = (x01 + x21)/2 
Interpolation Method 7: y7 = (xOO + x22)/2 
Interpolation Method 8: y8 = (x02 + x20)/2 
Interpolation Method 9: y9 = xOO 
45 Interpolation Method 10: y10 = x02 

Interpolation Method 11: y11 = x22 
Interpolation Method 12: y12 = x20 
Finally, the memory 15, the system controller 16 and the changing-over circuit 17 are described. 
The memory 15 stores data indicative of positions of a plurality of defective pixels contained in a solid im- 
50 aging device such as a CCD. The system controller 1 6 reads out the position data of the defective pixels from 
the memory 15 and controls the changing-over circuit 17 to produce the interpolation output when the position 
of the center pixel of the window of 3x3 taken in by the pixel taking-in circuit 1 0 is coincident with the position 
of the defective pixel stored in the memory 15 and to produce the center pixel x11 when not coincident. 

As described above, according to the defective pixel correction circuit of the present invention, the mag- 
55 nitudes of the boundaries between th peripheral or surrounding pixels are ordered to thereby identify the 
boundary pattern in the vicinity of the defective pixel so that the interpolation method suitable for the pattern 
is selected to perform interpolation. Accordingly, even when the complicated boundary is present in the vicinity 
of the defective pixel, the boundary can be identified to correct the defective pixel finely and exactly. 
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@ Verf^hren zum Betrieb eines digitalen Bildsystems einer Rontgendiagnostikeinrichtung 

(^) Die Erfindung betrifft ein Verfahren zum Betrieb eines 
digitalen Bildsystems (6) einer Rdntgendiagnostikeinrich- 
tung mit einer Rontgenelnheh (1, 2) zur Erzeugung von 
Rontgenbildern, einer Rdntgenbiidwandler-Fernsehkette (5 
bis 7) rur Erfassung der Rontgenbilder, die einen digitalen 
Bildwandler (5) mit matrixformig in Zeilen und Spalten 
angeordneten BUdpunkten und einen Monitor (7) zur Wie- 
dergabe der von dem digitalen Bild system (6) verarbeiteten 
Rontgenbilder aufweist. Zur Erkennung von defekten BUd- 
punkten wird wenigstens ein Kalibrierbild (10) durch eine 
HochpaS-Filterung (11), beispielsweise eine Median-Filte- 
rung (26 bis 29), in ein Hlterbild (12) gewandelt, das einer 
Defektbestimmung (13) zugefuhrt wird, so daB ein Defekt- 
bild (14) erhalten wird, das zur Korrektur (15) eines Original- 
bilds (16) verwendet wird. 
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1 

Beschreibung 

Die Erfindung betrifft ein Verfahren zum Betrieb ei- 
nes digitalen Bildsystems einer Rontgendiagnostikein- 
richtung mit einer R&ntgeneinheit zur Erzeugung von 5 
Rontgenbildern, einer Rontgenbildwandler-Fernsehket- 
te zur Erfassung der Rdntgenbilder, die einen digitalen 
Bildwandler mit matrixf ormig in Zeilen und Spalten an- 
geordneten BUdpunkten und einen Monitor zur Wieder- 
gabe der vom digitalen Bildsystem verarbeiteten Ront- 10 
genbilder aufweist. 

Zunehmend gewinnen beispielsweise die in der 
DE 43 00 828 Al beschriebenen digitalen bildgebenden 
Systeme in der Medizintechnik an Bedeutung, in denen 
sogenannte digitale Bildwandler die bisher eingesetzten 15 
analogen Bildwandler ersetzen. Derartig verwendete 
Bildwandler, beispielsweise die aus der DE 43 21 789 Al 
bekannten a-Si-Detektoren, lassen sich im allgemeinen 
noch nicht fehlerfrei produzieren. Vielmehr sind neben 
defekten Pixeln, bzw. BUdpunkten auch noch Gruppen- 20 
(Cluster-), Zeilen- und Spaltenausfalle zu verzeichnen. 
So werden beispielsweise Spaltenausfalle durch Unter- 
brechungen in den AdreDleitungen verursacht 

Um den AusschuB an Detektoren zu begrenzen, kdn- 
nen die defekten Bildpunkte korrigiert werden. Fiir ei- 25 
nen erfolgreiche Korrektur wird die Information bend- 
tigt, welche Pixel defekt und welche gut sind. 

Bisher wurden mehrere Dunkelbilder (Bilder ohne 
Strahlung) und mehrere Satze von Hellbildern, d. h. ob- 
jektlosen Bildern mit Strahlung, mit unterschiedlichen 30 
Intensitaten aufgenommen und durch eine arithmeti- 
sche Verknupfung ein Dunkelstrom-korrigiertes Bild 
erzeugt Aus diesen Bildern wurde dann versucht, die 
defekten Pixel zu finden. Dies konnte beispielsweise mit 
Histogrammen und bestimmten Schwellwerten im Hi- 35 
stogramm-Verlauf oder mit einer Summe von Histo- 
grammen fiir unterschiedliche Teilflachen des Bildes 
oder durch Prufung aller Bildpunkte des Bildes mit ge- 
wissen Bereichsschranken erreicht werden. 

AUe diese oder ahnliche Verfahren fuhren jedoch zu 40 
Problemen, wenn sich beispielsweise die elektronischen 
Gegebenheiten des Systems bei einem neuen Detektor 
oder^nderen Betriebsmodi des Detektors verandern, 
weil sie meist von einer Vielzahl von Parametern abhin- 
gen, die fur ein bestimmtes Detektor-Exemplar erst er- 45 
mittelt werden mussen. 

In einer Parallelanmeldung ist ein Verfahren zur De- 
fekterkennung bei einem a-Si-Panel beschrieben, bei 
der zur Erkennung von defekten Bildpunkten mehrfach 
auf zwei Bilder, von denen ein erstes Bild mit Rontgen- 50 
strahlung ohne Objekt und ein zweites Bild ohne Ront- 
genstrahlung erstellt ist, eine rekursive TiefpaB-Filte- 
rung nach-dem Prinzip der unscharf en Maske mit Tren- 
nung der Defekte nach Spalten und Zeilen angewandt 
wird. 55 

Die Erfindung geht von der Auf gabe aus, ein Verfah- 
ren der eingangs genannten Art zu schaffen, das auf 
einfache Weise eine pixelgenaue und sichere Erkennung 
defekter Bildpunkte ermdglicht 

Die Aufgabe wird erfindungsgem&B dadurch gel6st, eo 
daB zur Erkennung von defekten Bildpunkten wenig- 
stens ein Kalibrierbild durch eine Filterung in ein Filter- 
bild gewandelt wird, das einer Defektbestimmung zuge- 
fuhrt wird, so daB ein Defektbild erhalten wird, das zur 
Korrektur eines Originalbildes verwendet wird. 65 

ErfindungsgemaB kann die Filterung eine HochpaB- 
Filterung, beispielsweise eine zweidimensionale Me- 
dian-Filterung sein, bei der zuerst eindimensional in ei- 
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ner ersten Richtung und dann eindimensional in einer 
zweiten, zur ersten senkrechten Richtung gefiltert wer- 
den kann. 

Eine Bestimmung der Defekte kann erfindungsgemaB 
durch Ermittelung der Histogrammverteilung erfolgt. 

Es hat sich als vorteilhaft erwiesen, wenn folgende 
Verfahrensschritte angewandt werden: 

a) fuhre HochpaB-Filterung eines Kalibrierbildes 
durch, 

b) bestimme den Histogrammkern und 

c) ermittele Defekte der Bildpunkte auBerhalb des 
Histogrammkernes. 

Alternativ kann zur HochpaB-Filterung eine TiefpaB- 
Filterung des Kalibrierbildes durchgefuhrt und das der- 
art erhaltene Filterbild von dem Kalibrierbild subtra- 
hiert werden. 

Ein besonders einf aches und schnelles Verfahren er- 
halt man, wenn das Kalibrierbild einer Korrektur be- 
reits gefundener Defekte unterworfen wird, wenn das 
Korrekturbild durch eine TiefpaB-Filterung in ein Fil- 
terbild gewandelt wird, das von dem Kalibrierbild sub- 
trahiert wird, und wenn das derart erhaltene Differenz- 
bild zur Defektbestimmung einem Schwellenvergleich 
zugefuhrt wird, so daB ein Defektbild erhalten wird, das 
zur Korrektur des Kalibrierbild verwendet wird. 

Es hat sich als vorteilhaft erwiesen, wenn die Schwelle 
fur den Schwellenvergleich durch Bestimmung der Hi- 
stogrammverteilung des Differenzbildes ermittelt wird. 
ErfindungsgemaB kann die TiefpaB-Filterung nach- dem 
Prinzip der unscharf en Maske durchgefuhrt sein. ■ 

Es hat sich als vorteilhaft erwiesen, wenn folgende 
Verfahrensschritte angewandt werden: 

a) korrigiere im Originalbild bereits gefundene De- 
fekte, 

b) fuhre TiefpaB-Filterung des korrigierten Bildes 
durch, 

c) bilde die Differenz von dem Originalbild und 
dem Bild mit TiefpaB-Filterung, 

d) untersuche diese Differenz nach neuen Defekten 
und 

e) wiederhole den Vorgang, wenn neue Defekte 
gefunden wurden. 

Als Kalibrierbilder konnen erfindungsgemaB ein 
Hell- und ein Dunkelbild, ein Rausch- und ein Signalbild 
verwendet werden, wobei die Kalibrierbilder aus einer 
Serie von aufeinanderfolgender Kalibrierbilder beste- 
hen konnen. 

Eine genaue Ermittelung von weiteren Kalibriersi- 
gnalen kann erfolgen, wenn die Erkennung defekter 
Bildpunkte auch in einer Referenzzone durchgefuhrt 
wird. 

Als Rdntgenbildwandler mit defekten Bildpunkten 
kann erfindungsgemaB ein Bildwandler aus amorphem 
Silizium (aSi) oder ein Rontgenbildverstarker mit ange- 
koppeltem CCD- Bildwandler verwendet werden. 

Die Erfindung ist nachfolgend anhand von in der 
Zeichnung dargestellten Ausfuhrungsbeispielen naher 
erlautert Es zeigen: 

Fig. 1 eine Rontgendiagnostikeinrichtung nach dem 
Stand der Technik mit einem digitalen Bildsystem zur 
Anwendung des erfindungsgemaB en Verfahrens, 

Fig. 2 die allgemeinen Schritte des erfindungsgema- 
Ben Verfahrens, 

Fig. 3 eine erste Ausfiihrungsform des erfindungsge- 
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maBen Verfahrens und 

Fig. 4 eine alternative Ausfiihmngsform des erfin- 
dungsgemaBen Verfahrens. 

In der Fig. 1 ist eine bekannte Rontgendiagnostikein- 
richtung mit einer R6ntgenr6hre 1 dargestellt, die von 
einem Hochspannungsgenerator 2 betrieben wird. Die 
Rbntgenrohre 1 sendet ein Rontgenstrahlenbundel 3 
aus, das ein en Patienten 4 durchdringt und auf einen 
Rontgenbildwandler 5 entsprechend der Transparent 
des Patienten 4 geschwacht als RSntgenstrahlenbild 
fallt. Der Rontgenbildwandler 5 ist mit einem digitalen 
Bildsystem 6 und einen daran angeschlossenen Monitor 
7 2ur Wiedergabe des Rontgenstrahlenbildes verbun- 
den. Das digitale Bildsystem 6 kann in bekannter Weise 
eine Verarbeitungsschaltung, Wandler, Differenzstufen 
und Bildspeicher aufweisen. 

Bei einer derartigen RSntgendiagnostikeinrichtung 
kann der Rdntgenbildwandler 5 aus einer Photodioden- 
Matrix oder aus einem oder mehreren CCD-Bildsenso- 
ren bestehen, die mit einem RSntgenbildverstarker oder 
einer Szintillatorschicht aus einem gegenuber Rontgen- 
strahlen empfindlichen Leuchtstoff gekoppelt sincL Der 
Rontgenbildwandler 5 kann aus wasserstoffhalrigem 
amorphem Silizium (aSi : H) bestehen. 

Bei derartigen Bildwandlern konnen mehrere BOd- 
punkte der Bildpunkt-Matrix defekt sein, Es k6nnen ein- 
zelne Bildpunkte, mehrere Zeilen und/oder Spalten 
ubergreifende Gruppen von Bildpunkten (Cluster), gan- 
ze Zeilen oder Spalten defekt sein. 

Was sind nun^defekte Bildpunkte? Dazu muB sozusa- 
gen das "normale" Verhalten der Pixel bestimmt wer- 
den. UngleichmaBigkeiten in der Bestrahlung durch die 
R6ntgenr6hre 1 oder des verwendeten Szintilators bil- 
den sich normal erweise nur in niedrigen Ortsfrequenzen 
des Kalibrierbildes ab. Des weiteren sind die Rauschan- 
teile im Kalibrierbild unerwunscht und konnen unbe- 
rucksichtigt bleiben. Der erste Ansatz ware somit als 
Normal die mit einem nicht zu kleinem Kernel mit bei- 
spielsweise einer KernelgrSBen von 31 x31 besummte 
unscharfe Maske als Rechteckfilter des Kalibrierbildes 
zu definieren. Leider werden dabei die Defekte auf- 
grund ^aer Mittelung sozusagen auf ihre Umgebung ver- 
schmiert, wodurch auch Nachbarpixel mit beeinfluBt 
werden. 

Besser eignet sich ein Medianfilter fQr diesen FalL 
Allerdings weist dieser vor allem bei grfiBeren Kerneln 
eine recht hohe Rechenzeit auf. Diese laBt sich teilweise 
vermindem indem man den zweidimensionalen Filter in 
zwei eindimensionalen Filterungen aufteilt, die in Zei- 
len- und Spaltenrichtung getrennt nacheinander durch- 
gefiihrt werden. Der Unterschied gegenuber dem ech- 
len zweidimensionalen Filter ist gering und fOr diesen 
Zweck vemachlassigbar. 

Anhand der Fig. 2 wird das allgemeine erfindungsge- 
rnaBe Verfahren dargestellt Ein Kalibrierbild 10 wird 
einer HocbpaB-Filterung 11 unterworfen. Das erhaltene 
rilterbild 12 wird einer Defektbestimmung 13 zuge- 
fiihrt, als deren Ergebnis ein Defektbild 14 vorliegt, das 
samtliche defekte Bildpunkte kennzeichnet Mit diesem 
Defektbild 14 kann eine Konrektur 15 eines driginalbil- 
des 16 mit Rontgenstrahlung erfolgen, so daB man ein 
korrigiertes Rontgenbild 17 erhalt. Die Korrektur 15 
kann beispielsweise dadurch erreicht werden, daB die 
defekten Bildpunkte durch eine lineare Interpolation 
benachbarter Bildpunkte ersetzt wird. 

In Fig. 3 ist nun ein erstes erfindungsgemafies Verfah- 
ren wiedergegeben. Es werden vier verschiedene Bilder 
aJs KalibrierbOder 10 verwendet, ein Rausch-Dunkelbild 
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20, ein Signal -Dunkelbild 21, ein Signal-Hellbiid 22 und 
ein Rausch-Hellbild 23. Dabei bedeutet Dunkelbild ein 
BDd bzw. eine Bildserie ohne jegliche Belichtung und ein 
Hellbild ein Bild bzw. eine Bildserie mit gleichformiger 
Belichtung ohne Untersuchungsobjekt Als Rausch-Bild 
wird die Streuung oder Varianz a2 jeweils fur jedes 
Pixel einer Bildserie und als Signal ebenfalls jeweils fur 
jedes Pixel der Mittelwert einer Serie von Bildern ge- 
wertet Auch ein Signal eines Dunkelbildes der Refe- 
renzzone 24 wird entsprechend erf aBt 

Die Dunkel-Referenzzone des verwendeten a-Si-De- 
tektors am linken Bildrand, dem Zeilenanfang, bedarf 
einer besonderen Behandlung. Diese ist licht- und strah- 
lungsdicht abgedeckt, so daB auch im Hellbild kein zu- 
satzliches Signal entsteht Sie dient dazu, das Zeilenrau- 
schen zu vermindem und kleine Transienten im Offset 
zu korrigiereiL Auch hier mussen Defekte erkannt wer- 
den, damit diese nicht die Zeilenkorrektur verfalschen. 
Als Kriterium beziiglich des Signals muB hier der (Dun- 
kelstrom-)Offset alleine dienen. Die Defektschranken 
konnen hier sogar etwas enger gesteckt werden, da eini- 
ge zusatzliche Defekte nicht die Bildinformation schma- 
lern, auf der anderen Seite aber wQrden hier nicht er- 
kannte Defekte die Minderung des Zeilenrauschens zu- 
25 nichte machen. Zeilen mit zuwenig guten Pixeln werden 
als ganzes als defekt erkl art- 
Die beiden Rauschbilder 20 und 23, das Signalbild der 
Referenzzone 24 sowie die Differenz 25 der Signalbilder 
21 und 22 werden einer zweidimensionalen HochpaB- 
Filterung, einer Median-Filterung 26 bis 29 unterwor- 
fen. Bei einer derartigen Median-Filterung wird im Ge- 
gensatz zum einfachen Rechteckfilter in der Umgebung 
eines Bildpunktes nicht der Mittelwert gebildet, sondern 
wird der, nachdem alle in der Umgebung gefundenen 
Werte sortiert wurden, in der Mitte liegende Wert als 
Filterwert genommen. Einzelne, weit vom Mittelwert 
entfernte Werte, die wahrscheinlich Defekte sind, haben 
so praktisch keinen EinfluB auf den Filterwert. So wie 
beschrieben hat der Filter einen TlefpaB-Charakter. 
Den HochpaB-Charakter erhalt man durch einfache 
Subtraktion des Ausgangswertes vom Originalwert. 

Aus den Filterwerten wird die Histogranunverteilung 
30 bis 33 ermittelt und anschlieBend aus diesem Histo- 
gramm die Defektbestimmung 34 bis 37 durchgefiihrt. 
Dies erfolgt derart, daB jeweils der Kern der Histo- 
grammverteilung 30 bis 33 ermittelt wird, das heLBt, der 
Bereich, der sich um das Maximum der gausformigen 
Kurve erstreckt. Durch einfache Erweiterung dieses Be- 
reiches um einen Faktor werden die defekten Bildpunk- 
te ermittelt, die auflerhalb dieses Bereiches liegen. 

Die Ausgangssignale der Defektbestimmungen 34 bis 
37 werden einer Oder-Verknupfung 38 zugefuhrt, an 
deren Ausgang das Defektbild 14 liegl 

Anhand der Fig. 4 wird nun eine Alternativlosung des 
erfindungsgemaBen Verfahrens beschrieben, die sich 
insbesondere durch geringen Rechenaufwand auszeich- 
net. Die Median-Filterung wird hier durch eine Iteration 
ersetzt, die den Rechteckfilter benutzt. Ansonsten blei- 
ben die Verarbeitungsschritte gleich. Das Kalibrierbild 
10, das jeweilige Signal- oder Rauschbild 20 bis 24, wird 
einer Defektkorrektur 40 zugefuhrt, in der bereits er- 
kannte defekte Bildpunkte beispielsweise mit einer li- 
nearen Interpolation korrigiert werden. Mach einer 
TiefpaB-Filterung 41 mit einem Rechteckfilter wird das 
Filterbild 42 einer Subtraktion 43 mit dem Kalibrierbild 
10 unterworfen. Durch diese TiefpaB-Filterung 41 und 
anschlieBende Subtraktion 43 erhalt man ein gefdtertes 
Signal, ahnlich wie durch die HochpaB-Filterung 1 1. Aus 
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dem Subtraktionsbild 44 werden uber einen Schwellen- 
vergleich 48 neue Defekte gefunden, die in das Defekt- 
bild 14 eingetragen werden. Dieser Vorgang wird mehr- 
fach wiederholt 

Die TiefpaB-Filterung 41 durch ein Rechteckfilter mit 5 
einem nicht zu kieinem Kernel durchgefuhrt, wobei das 
Rechteckfilter um einen Bildpunkt einen arithmetischen 
Mittelwert bildet 

Zur Einstellung der Schwelle wird aus dem Subtrak- 
tionsbild 44 das Histogramm 45 ermittelt, aus dem der 10 
Kernbereich 46 bestimmt wird Die Bereichserweite- 
rung 47 um einen Faktor bestimmt dann die Hone der 
SchweUe fur den Schwellenvergleich 48 der defekten 
Bildpunkte. 

Diese Iteration zur Bestimmung der defekten Bild- 15 
punkte weist gegenuber der Median-Filterung geringe- 
re Rechenzeiten auf. Der Vorteil der Median-Filterung 
jedoch ist, daB man mit einem Schritt das gewunschte 
Def ektbild erhalt 

Die Rechenzeiten der Median-Filterung lassen sich 20 
durch folgende MaBnahme beschleunigen. 

Bei der Orginal-Median-Filterung muB der Kern um 
jedes Pixel sortiert werden, d. h. im Falle eines 30*30 
Kerns in einem 1000*1000 Bild muB 106 rnal ein Feld von 
900 Werten sortiert werden, so daB sich zusammen 25 
900*106 Sortierwerte ergeben. 

Wird dagegen die Median-Filterung nach x und y auf- 
geteilt und zweimal hintereinander eine Median-Filte- 
rung gebildet, indem einmal mit einem Kern von 30*1 
und dieses Ergebnis dann mit einem Kern von 1 *30 gefil- 30 
tert wird. Dann sind nur 2*106 Felder von 30 Werten zu 
sortieren, so daB sich 60*106 Sortierwerte ergeben, so 
daB sich 15 mal weniger Sortierwerte als vorher erge- 
ben. 

Will man trotzdem, wegen der schnelleren Berechen- 35 
barkeit, mit der unscharfen Maske arbeiten, muB man 
den EinfluB der Defekte auf ihre Nachbarschaft vermin- 
dern. Dies wird erfindungsgemiB durch folgende Itera- 
tion gel6st: 

40 

a) korrigiere im Kalibrierbild bereits gefundene 
Iiefekte, beispielsweise mit linearer Interpolation, 
Ze z. B. mit einfachen Schwellen gefunden wurden, 
die ein Signal nahe der Bereichsgrenzen des A/D- 
Wandlers haben, da die Kalibrierbiider normaler- 45 
weise weit innerhalb dieser Grenzen aufgenommen 
werden, 

b) fQhre eine TiefpaB-Filterung des korrigierten 
Bildes durch, 

c) bilde Differenz vom ursprunglichen Kalibrierbild 50 
und TiefpaB-gefilterten, korrigierten Bild, 

d) bestimme aus dem Histogramm der Differenz 
Schwellen fiir Defekte, wobei im Histogramm be- 
reits gefundene Defekte unberucksichtigt bleiben, 

e) finde neue Defekte und fuge sie den bereits ge- 55 
fundenen hinzu und 

f) wiederhole Iteration bis kaum noch neue Defekte 
gefunden werden. 

Bezuglich der Schwellwertftndung kann beispielswei- 60 
se angenommen werden, daB der Kernbereich des Hi- 
stogramms in einer Normalverteilung, der 50% aller 
Pixel enth&lt, 0,675 Mal der Standardabweichung ent- 
spricht Bei einer Erweiterung dieses Bereiches um bei- 
spielsweise einen Faktor 7 sind bei einer Normalvertei- 65 
lung nur noch 3*10-5 aller Werte auBerhalb zu finden. 
Pixel, die trotzdem noch auBerhalb zu finden sind, sind 
mit hoher Wahrscheinlichkeit defekt 
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Das Histogramm nahert sich durch die Iteration der 
Normalverteilung an, da immer mehr bereits gefundene 
Defekte ausgeschlossen und im Histogramm nicht mehr 
berucksichtigt werden. 

Der Vorteil dieser Methode liegt in der automati- 
schen Anpassung der Schwellen an die im Kalibrierbild 
vorhandenen normalen Rauschwerte der PixeL 

Die Verfahren werden erfindungsgem&B bei einem 
a-Si-Detektor fur die Kalibrierung verwendet Ein der- 
artiger Detektor kann beispielsweise eine Matrix von 
1024 Zeilen und 1080 Spalten aufweisen. Die Referenz- 
zone kann aus den Bildpunkten Nr. 0 bis 79 bestehen, 
wobei beispielsweise Nr. 1 bis 64 davon verwendet wer- 
den und mindestens 16 nicht defekt sein sollten. Der 
A/D-Wandler kann eine Tiefe von 12 Bit aufweisen. 

Zur Kalibrierung werden immer eine Serie Dunkel- 
und eine Serie Hellbilder aufgenommea Folgende 
Schritte werden angewendet: 

— Alle Bilder werden auf Einhaltung der Grenzen 
des A/D-Wandlers untersucht Da bei der Aufnah- 
me der Kalibrierbiider sichergestellt sein muB, daB 
nirgends eine Oberschreitung des Digitalisierungs- 
bereichs des A/D-Wandlers auftritt, sind Pbcel, die 
trotzdem ein Signal nahe der Grenzen des A/D- 
Wandlers aufweisen, sicher defekt Diese Defekte 
sind ein guter Ausgangspunkt fiir die oben be- 
schriebenen Iteration. 

— Die Bilder jeder Serie werden gemittelt, sowie 
das Rauschen in den Pbcel berechnet wird. Die vier 
resultierenden Bilder werden der oben beschriebe- 
nen Defekterkennung zugefuhrt Bei kurzen Kali- 
brierserien kann man das Rauschen Zeilen- und 
Spaltenweise mitteln. 

— Nur in der Referenzzone stutzt sich die Defekt- 
erkennung nur auf die Offsetwerte, da durch die 
Abdeckung kein anderweitiges Signal erzeugt wer- 
den kann. Defekte werden mittels der oben be- 
schriebenen Median-Filterung bzw. Iteration er- 
kannt Zeilen mit zu vielen Defekten in der Refe- 
renzzone werden ausgesondert Einige irrtumlich 
als defekt erkannte Pbcel sind hier nicht so gravie- 
rend Anders im eigentlichen Bild, wo mdglichst viel 
Information erhalten bleiben sollte. 

Alle vier Bilder konnen mit dem oben angefuhrten 
Verfahren auf Defekte untersucht werden. Allerdings 
wird man beim Signal besser nur die Differenz von Hell- 
und Dunkelbild untersuchen, da der (Dunkelstrom-)Of f- 
set allein kein gutes Defektkriterium sein muB, da er 
spater bei der Aufnahme leicht korrigiert werden kann. 

Als wichtigster Schritt bei der Defekterkennung er- 
wies sich die Untersuchung der Rauschwerte, da in Bild- 
folgen stark rauschende oder blinkende Pbcel oder gar 
flackernde Zeilen und Spalten enorm den Bildeindruck 
storen. Beim Signal fallen diese nicht unbedingt auf, da 
sich durchaus ein normaler Mittelwert einstellen kann. 
Auch eine Defekterkennung, die sich im wesentlichen 
nur auf die Rauschwerte stutzte, zeigte sehr gute Ergeb- 
nisse. 

Bei kurzen Serien kann die Statistik uber die Rausch- 
werte recht schlecht sein. Hier bietet sich dann an, uber 
bestimmte Pixelgruppen bis hin zu Zeilen und Spalten 
zu mittem. Allerdings k6nnen so keine einzelnen Defek- 
te mehr erkannt werden. Andererseits kann nach Mitte- 
lung uber Zeilen/Spalten eindimensional weitergearbei- 
tet werden, so daB der Medianfilter wieder anwendbar 
ist 
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Gegenilber fruheren Verfahren, die sich hauptsach- 
Hch auf die Signalwerte stfltzten, wurden mit dem erfin- 
dungsgemaBen Verfahren wesentlich weniger Bild- 
punkte unnfitigerweise als defekt erkannt, wobei sich 
aber der Bildeindruck wesentlich verbesserte. Letzteres 5 
ist wohl der zusatzlichen Rauschuntersuchung zuzu- 
schreiben. 

Der Detektor wird in mehreren, unterschiediich emp- 
findlichen Modi betrieben, was auch EinfluB auf das 
Rauschen der Pixel haL Trotzdem zeigte sich das Ver- 50 
fahren unempfindlich dagegen, da alle Schwellen auto 
matisch berechnet werden, bzw. relativ zu berechneten 
Werten eingestellt werden* 

Alle beschriebenen Verfahrensschritte k6nnen anstel- 
le durch eine das Bildsystem 6 steuemde Software auch 15 
als Hardware-L6sung durchgefuhrt werden. 

Die Erkennung und Korrektur defekter Bildpunkte 
kann effizient in Hardware realisiert oder von einem 
Multi-Prozessorsystem ausgefuhrt werden, wobei hier 
die Einzel-Prozessoren optimal einsetzbar sind 20 

Patentanspriiche 

1. Verfahren zum Betrieb eines digitalen Bildsy- 
stems (6) einer Rdntgendiagnostikeinrichtung mit 25 
einer Rfcntgeneinheit (1, 2) zur Erzeugung von 
Rdntgenbildern, einer R6ntgenbildwandler-Fern- 
sehkette (5 bis 7) zur Erfassung der Rontgenbilder, 
die einen digitalen Bildwandler (5) mit matrixfdr- 
mig in Zeilen und Spaiten angeordneten Bildpunk- 30 
ten und einen Monitor (7) zur Wiedergabe der von 
dem digitalen Bildsystem (6) verarbeiteten Ront- 
genbilder aufweist, dadurch gekennzeichnet, daB 
zur Erkennung von defekten Bildpunkten wenig- 
stens ein Kalibrierbild (10) durch eine Filterung (11) 35 
in ein Filterbild (12) gewandelt wird, das einer De- 
fektbestinunung (13) zugefuhrt wird, so daB ein De- 
fektbild (14) erhalten wird, das zur Korrektur (15) 
eines Originalbild (16) verwendet wird 

2. Verfahren nach Anspruch 1, dadurch gekenn- 40 
zejchnet, daB zur Filterung (11) eine HochpaB-Fil- 
terung mit einem Median-Filter (26 bis 29) durchge- 
fuhrt wird 

3. Verfahren nach Anspruch 1 oder 2, dadurch ge- 
kennzeichnet, daB eine zweidimensionale Median- 45 
Filterung (26 bis 29) zuerst eindimensional in einer 
ersten Richtung und dann eindimensional in einer 
zweiten, zur ersten senkrechten Richtung durchge- 
fuhrt wird. 

4. Verfahren nach einem der AnsprOche 1 bis 3, 50 
dadurch gekennzeichnet, daB die Defektbestirn- 
mung (13) durch Ermittlung der Histogrammvenei- 
lung (30 bis 33) erf olgt. 

5. Verfahren nach einem der Anspruche 1 bis 4, 
gekennzeichnet durch folgende Verfahrensschritte: 55 

a) fuhre HochpaB-Filterung eines Kalibrierbil- 
des durch, 

b) bestirnme den Histogrammkern und 

c) ermittele Defekte der Bildpunkte aufierhalb 
des Histogrammkern es. 60 

6. Verfahren nach Anspruch 1, dadurch gekenn- 
zeichnet, daB zur HochpaB-Filterung (1 1) eine Tief- 
paB- Filterung (41) des Kalibrierbildes (10) durchge- 
fuhrt und das derart erhaltene Filterbild (42) von 
dem Kalibrierbild (10) subtrahiert wird 55 

7. Verfahren nach Anspruch 1 oder 6, dadurch ge- 
kennzeichnet, daB das Kalibrierbild (10) einer Kor- 
rektur (40) bereits gefundener Defekte unterwor- 
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fen wird, daB das Korrekturbild durch eine TiefpaB- 
Filterung (41) in ein Filterbild (42) gewandelt wird, 
das von dem Kalibrierbild (10) subtrahien wird, 
und daB das derart erhaltene Differenzbild (44) zur 
Defektbestimmung einem Schwellenvergleich (48) 
zugefiihrt wird, so daB ein Defektbild (14) erhalten 
wird, das zur Korrektur (40) des Kalibrierbild (10) 
verwendet wird 

8. Verfahren nach Anspruch 7, dadurch gekenn- 
zeichnet, daB die Schweile fur den Schwellenver- 
gleich (48) durch Bestimmung der Histogram rnver- 
teilung (45) des Differenzbild es (44) ermittelt wird 

9. Verfahren nach einem der Anspruche 6 bis 8; 
dadurch gekennzeichnet, daB die HefpaB-Filterung 
(41) nach dem Prinzip der unscharfen Maske durch- 
gefuhrt wird 

10. Verfahren nach einem der Ansprilche 6 bis 9, 
gekennzeichnet durch folgende Verfahrensschritte: 

a) korrigiere im Originalbild bereits gefundene 
Defekte, 

b) fuhre TiefpaB-Filterung des korrigierten 
Bildes durch, 

c) bilde die Differenz von dem Originalbild und 
dem Bild mit TiefpaB-Filterung 

d) untersuche diese Differenz nach neuen De- 
fekten und 

e) wiederhole den Vorgang, wenn neue Defek- 
te gefunden wurden. 

11. Verfahren nach einem der Anspruche 1 bis 10, 
dadurch gekennzeichnet, dafl als Kalibrierbild er 
(10) ein Hell- (22, 23) und ein Dunkelbild (20, 21) 
verwendet werden. 

12. Verfahren nach einem der Anspruche 1 bis 11, 
dadurch gekennzeichnet, daB als Kalibrierbilder 
(10) ein Rausch- (20, 23) und ein Signalbild (21, 22) 
verwendet werden. 

13. Verfahren nach einem der Anspruche I bis 12, 
dadurch gekennzeichnet, daB eine Serie von aufein- 
anderfolgenden Kalibrierbildern (10) verwendet 
wird 

14. Verfahren nach einem der Anspruche 1 bis 13, 
dadurch gekennzeichnet, daB die Erkennung defek- 
ter Bildpunkte auch in einer Referenzzone (25) 
durchgefuhrt wird 

15. Verfahren nach einem der Anspruche 1 bis 14, 
dadurch gekennzeichnet, daB der R6ntgenbild- 
wandler (5) mit defekten Bildpunkten ein Bildwand- 
ler aus amorphem Silizium ist. 

16. Verfahren nach einem der Ansprilche 1 bis 14, 
dadurch gekennzeichnet, daB der Rdntgenbild- 
wandler (5) mit defekten Bildpunkten ein Rontgen- 
bildverstarker mit angekoppeltem CCD-Bildwand- 
ler ist. 
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Claims 

1. Ad fective pixel correction circuit comprising: 

a pixel taking-in circuit (10) for taking in signals produced from peripheral pixels of a defective pixel 
in an imaging device; 

a boundary detection circuit (11) for detecting magnitudes of boundaries between adjacent pixels 
of the taken-in peripheral pixels; 

a boundary ordering circuit (12) using outputs of said boundary detection circuit for ordering the 
boundaries; 

an interpolation method determining circuit (13) for determining an interpolation method on the ba- 
sis of outputs of said boundary ordering circuit, and 

an interpolation circuit (14) for producing an interpolation signal by said determined interpolation 
method. 

2. A defective pixel correction circuit comprising: 

a pixel taking-in circuit (10) for taking in signals produced from peripheral pixels of a defective pixel 
in an imaging device; 

a boundary detection circuit (11) for detecting magnitudes of boundaries between adjacent pixels 
of the taken-in peripheral pixels; 

a boundary ordering circuit (12) using outputs of said boundary detection circuit for ordering the 
boundaries; 

an interpolation method determining circuit (13) for determining an interpolation method on the ba- 
sis of outputs of said boundary ordering circuit; 

an interpolation circuit (14) for producing an interpolation signal by said determined interpolation 
method; 

a memory (15) for storing a position of said defective pixel; and " 
replacement control means (17) for judging whether data is stored in said memory or not and re- 
placing an image of said pixel taking-in circuit with an image of said interpolation circuit when data is stored 
in said memory. 

3. A defective pixel correction circuit according to Claim 1, wherein said boundary detection circuit is to cal- 
culate the magnitudes of the boundaries between the adjacent pixels from eight peripheral pixels in said 
pixel taking-in circuit and said boundary ordering circuit is to effect determination on the basis of the mag- 
nitudes of the boundaries between the pixels. 

4. A defective pixel correction circuit according to Claim 2, wherein said boundary detection circuit is to cal- 
culate the magnitudes of the boundaries between the adjacent pixels from eight peripheral pixels in said 
pixel taking-in circuit and said boundary ordering circuit is to effect determination on the basis of the mag- 
nitudes of the boundaries between the pixels. 

5. A defective pixel correction circuit according to Claim 3 or 4, wherein said boundary ordering circuit in- 
cludes a rule forjudging that the magnitude of the boundary in the vertical direction is larger than that in 
the horizontal direction when the magnitudes of the boundaries in the vertical and horizontal directions 
are equal to each other. 

6. A defective pixel correction circuit according to Claim 3 or 4, wherein said boundary ordering circuit in- 
cludes a rule for judging that the magnitude of the boundary having a smaller suffix is larger than that 
having a larger suffix when the magnitudes of the boundaries in the vertical direction or the horizontal 
direction are equal to each other. 

7. A defective pixel correction circuit according to Claim 3 or 4, wherein said boundary ordering circuit in- 
cludes a rule forjudging that the magnitude of the boundary in the vertical direction is larger than that in 
the horizontal direction when the magnitude of the boundary in the vertical direction is equal to that in 
the horizontal direction and judging that the magnitude of the boundary having a smaller suffix is larger 
than that having a larger suff be when the magnitudes of the boundaries in the vertical direction or the hor- 
izontal direction are equal to each other. 
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